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The relaxation in complex systems is in general nonexponential. After an initial rapid decay 
the system relaxes slowly following a long time tail. In the present paper a sandpile modelation of 
the relaxation in complex systems is analysed. Complexity is introduced by a process of avalanches 
in the Bethe lattice and a feedback mechanism which leads to slower decay with increasing time. 
In this way, some features of relaxation in complex systems: long time tails relaxation, aging, and 
fractal distribution of characteristic times, are obtained by simple computer simulations. 
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About 150 years ago Weber and Gauss carried out a 
simple experiment demonstrating that relaxation in com- 
plex systems is not exponential. Investigating the con- 
traction of a silk thread they found that it does not con- 
tract so quickly but it relaxes slowly following a power 
law t~ a . This behavior is not particular of mechanical 
relaxation, it has been also observed in experiments of 
magnetic relaxation in spin glasses and high-Tc supercon- 
ductors, transient current measurements in amorphous 
semiconductors, dielectric relaxation, and more Q. This 
phenomenon of slow relaxation following a power law is 
called a long time tail and it seems to be characteristic 
of the relaxation in complex systems. 

Associated to the slow relaxation dynamics aging is 
observed, which means that the properties of a system 
depends on its age ||. For instance, consider a glass 
quenched at time t = below its glass transition tem- 
perature under an external stress and at time t = t w the 
stress is released. If the system were near equilibrium 
then its response measured at certain time t > t w will 
be a function of the difference t — t w and, therefore, nor- 
malized responses taken at different initial times t w will 
collapse in a single curve. However, in the glass mate- 
rial it is observed that the relaxation becomes slower for 
larger t w . The aging phenomena is a consequence of the 
nonequilibrium dynamics and may be also considered a 
characteristic of the dynamics of complex systems out 
from equilibrium. 

Sandpiles seem to be simplest systems which also lead 
to complex behavior.. They have been taken as the 
paradigm of self-organization since the introduction of 
these ideas by Bak, Tang, and Wiesenfeld ||. The evo- 
lution of a sandpile under an external perturbation has 
been extensively studied experimentally theoreti- 



cally [^,0, and by computer simulations M. In general 
these works analyze the sandpile dynamics close to the 
critical state, i.e. close to the critical angle. 

Slow relaxation dynamics have been observed in the 
process of compactation of sand heaps under vibrations 
UPl • On the other hand, Jaeger et al have measured the 
temporal decay of the slope angle of sandpiles under vi- 
brations with different intensities B . They observed that 
for high vibration intensities the pile angle relaxes loga- 
rithmically slowly to a steady state but it relaxes slower 
for low intensities. In order to explain the slow dynamics 
they proposed a simple model, very similar to flux creep 
models for magnetic relaxation in superconductors. In 
this way they explain the logarithmic decay observed for 
high vibration intensities but they could not account for 
the behavior at low intensities. The model proposed by 
Jaeger et al Q] fails at low vibration intensities because it 
does not take into account cooperative phenomena which 
leads to formation of avalanches. The avalanche dynam- 
ics should be the main contribution to the relaxation of 
the pile angle at low vibration intensities, close to the 
static limit. 

In the present paper the angle relaxation of a sand- 
pile under vibrations is investigated. Avalanches are 
modelled by a mean-field approximation represented in a 
Bethe lattice. Typical characteristics of complex systems: 
long time tails relaxation, aging and fractal distributions 
of time constants are obtained by simple numerical sim- 
ulations. The main goal of this work if to study the 
sandpile dynamics far from equilibrium. 

In a static sandpile sand grains are in equilibrium un- 
less the angle of the pile overcomes its critical value de- 
termined by the static friction coefficient. However, the 
grains in a pile under vibration are vibrating around cer- 
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tain quasiequilibrium positions but from time to time the 
amplitude of its vibrations becomes large. These large 
amplitude motions carry as a consequence a convective 
motion in the bulk and the relaxation of the pile angle. 

At low intensity vibrations large amplitude grain mo- 
tions are rare events but from time to time a grain may 
jump from its quasiequilibrium position. If a surface 
grain jumps it will fall through the slope of the pile un- 
til it collides with another sand grains down the slope. 
After this collision the initial grain may be trapped with 
those grains or some of them may fall through the slope. 
If the last possibility happens then each grain "surviv- 
ing" the collision will fall through the slope as the initial 
grain did. This image of an avalanche as an initial object 
that consecutively drags another resembles a branching 
process for which the Bethe lattice representation seems 
to be natural |?]]. 

However, representing the avalanche dynamics in a 
Bethe lattice is a mean field approach to the prob- 
lem, it neglects correlations between branches. Notwith- 
standing, if the pile angle is below is critical value the 
avalanches will be rare events and, in case of occurrence, 
will be very sparse. Thus, such an approximation will be 
acceptable for analyzing the long time relaxation of the 
pile angle below its critical value, which is the subject of 
this paper. Away from the critical point the mean field 
approximation works quite well. 

Let us represent the avalanche as a cascade in the 
Bethe lattice as follows. Firstly, we start with a sin- 
gle node, which could represent in this case a grain. In a 
further step F will emerge with probability p(8), depend- 
ing on the pile angle 9. This operation of generation of 
F identical particles starting from one is repeated in the 
next step to each node of the new group, and so on. If 
the percolation process overcomes a given length (that of 
the border of the pile) those nodes beyond the limit con- 
stitute the avalanche. If it does not, there would be no 
avalanche since the cascade was stopped before reaching 
the base of the pile (frustrated avalanche). By avalanche 
size we take the number of nodes, of the corresponding 
Bethe lattice, in the last step. 

After an avalanche the sandpile autorganizes itself with 
the new number of grains (i.e. a new slope is calculated 
with the remaining grains). The occurrence of avalanches 
will carry as consequence a decrease in the number of 
grains in the pile and, therefore, a decrease of p(9). As 
more avalanches take place less probable is the occurrence 
of a new avalanche. Thus to characterize this feedback 
mechanism a relation between p and the number of grains 
in the pile N is needed. 

The drag probability p(9) is a function of the slope. Its 
value is determined by the competition of two contrary 
forces: the gravity which conspires against the stability 
of the slope and the friction which favours the slope sta- 
bility. Since the slope forms an angle 9 with the horizon- 
tal plane the component of the gravity force in the slope 
direction will be larger with the increase of this angle, 
varying from zero to a maximum value when 9 goes from 



zero to 7r/2. Therefore it is plausible to assume that the 
contribution of the gravity and, therefore, the tendency 
of falling down the slope is proportional to sin^. On 
the contrary, the resistance to this tendency given by the 
static friction decreases with decreasing the pile angle, 
varying from a maximum value to zero when 9 goes from 
zero to 7r/2. Hence, it is also plausible to assume that 
the resistance to the falling down is proportional to cos 9. 

Thus the slope dependence of p should be given 
through the ratio of both tendencies sin 9/ cos9 = tan#. 
Based on this hypothesis we propose the exponential re- 
lation 

p(9) = exp(- J 4/tan6») , (1) 

where A is a parameter determined by the gravitational 
field, the friction, and vibration intensity. Notices that 
p(0) = and p(tt/2) = 1. 

On the other hand, as it is well known the number of 
grains in a pile with slope angle 9 is given by 

N = Bx 3 tan 9 (2) 

where B is a geometrical factor and x = D/d is the ra- 
tion between a characteristic size of the pile base D and 
sand grain d. Combining equations [l] and ^| it is obtained 

p(N) = cxp{-cx 3 /N) . (3) 

where c = AB. Thus equation |l] relates the dragging 
probability with the number of grains in the pile. 

The parameter c should not depend on the size of the 
system. It must be a function of parameters describing 
the vibrations, amplitude a and frequency uj, and of the 
gravitational field acceleration g. The only nondimen- 
tional combination of this magnitudes is given by the 
ratio between the vibration acceleration aoj 2 and g and, 
therefore, c = c(au 2 /g). This conclusion obtained from 
dimensional analysis is corroborated by experiments on 
sandpiles under vibrations fll| which shown that the ra- 
tio au 2 / g is the relevant parameter. 

For a given value of c and x there is a critical number 
of grains in the pile N c . This value can be found recall- 
ing that in the Bethe lattice the critical value of p for 
percolation exists is l/F, resulting 

iV c = cx 3 /lnF. (4) 

In a static sandpile (i.e. no vibrations) this value corre- 
sponds with the number of grains in the pile at the angle 
of repose . 

Another magnitude of interest is the penetration 
length, the number of steps n{9) in the Bethe lattice for 
an avalanche take place. This magnitude must be pro- 
portional to the length of the pile slope and, therefore, 
should be given by the expression 

n(9) = x/ cos 6 (5) 

Notice that the possible existence of a geometrical fac- 
tor in equation || may be absorbed in x, redefining B in 
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equation Moreover, the relation between n and N may 
be easily obtained using equation ^[ 

0.0 - 



Equations || and || relates the parameters of the pile 
with those of its Bethe lattice representation. They were 
obtained here in a different way than in . Other depen- 
dencies between the dragging probability and the num- 
ber of grains in the sandpile may be proposed. Notwith- 
standing, as it is discussed latter, the precise form of this 
functional dependence is not too relevant. 



The numerical experiment of relaxation is performed 
as follows. We start with a certain number of N. In a 
first step, we test if an avalanche takes place using the 
Bethe lattice representation. If it does then we recalcu- 
late the value of p(N) by simply resting the size of the 
avalanche to N and using equation 0. Then this step is 
repeated again and again. 
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FIG. 2. Normalized relaxation <j>(t) = N(t)/N(t = 0) of 
the number of grains in the pile for different values of c. The 
initial condition is N(t = 0) = N c in all cases. Insets show 
the values of a obtained from the fit to the power tail t~ a . 



If avalanches are considered as instantaneous the num- 
ber of steps is a measure of time. This approxima- 
tion is valid for low vibrations intensities. In this case 
grain jumps which trigger avalanches are rare events and, 
therefore, the time between to successive grain jumps will 
be much larger than the duration of avalanches. 
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FIG. 1. Normalized relaxation 4>{t) = N(t)/N(t = 0) of 
the number of grains in the pile for different values of x. The 
initial condition is N(t = 0) == N c in all cases. Insets show 
the values of a obtained from the fit to the power tail t~ a . 



The normalized decay of N in time for different values 
of x and c is shown in figures Qand|. In all cases we took 
as initialcondition N(t = 0) = N c , the number of grains 
at the angle of repose in the corresponding static limit, 
and we have averaged over 10000 realizations. As it can 
be seen the relaxation follows a fast decay at short times 
but it is slower than an exponential at long times. The 
log-log plot reveals a straight line at long times charac- 
teristic of a long time tail t~ a . Insets show the values of 
a obtained from the fit to this power tail, having relative 
errors of the order of 0.01. 



The characteristic exponent a decreases towards zero 
with increasing x, i.e. the relaxation is slower for larger 
sandpiles. It is expected that it is strictly zero for x — > oo, 
which means that an infinite sandpile will be in equi- 
librium. This limit may be the equivalent of the ferro- 
magnetic state of magnetic materials where domain are 
macroscopic structures having finite magnetization. On 
the other hand, with increasing c, a decreases towards 
zero and since the relaxation should be slower for smaller 
vibrations intensities we may conclude that c must in- 
crease with decreasing the vibration intensity. Hence, c 
is a monotically decreasing function of auj 2 /g. 
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FIG. 3. Aging. Normalized 

relaxation function <j>(t + t w ,t w ) = N(t + t w )/N(t w ) taken 
at three different initial times t w . 
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FIG. 4. Distribution of time between avalanches. The 
dashed line is a power tail with exponent — 1 — a, a being 
the exponent of the long time tail shown in figure 1. 



The existence of aging in our model is illustrated in 
figure |^. We have plotted the normalized relaxation at 
different ages (steps) of the system(simulation) by tak- 
ing as initial time t — t w . As it can be see the re- 
laxation becomes slower with increasing the age of the 
system, in agreement with experiments in structural and 
spin glasses Q . Thus, since the angle relaxation becomes 
slower with the age of the pile then it will never reach an 
equilibrium angle and, therefore, properties like transla- 
tional invariance and the fluctuation dissipation theorem 
do not hold §. 

Associated with these slow relaxation dynamics and 
aging phenomena we expect to observe a wide distri- 
bution of time between avalanches. As time increases 
the time between two consecutive avalanches At becomes 
larger and larger, because after an avalanche the occur- 
rence of a new avalanche becomes smaller. Therefore, 
it is expected that the mean time between avalanches 
diverges as t — > oo. Hence, the distribution of time be- 
tween avalanches n(At) should satisfies the asymptotic 
behavior for large At 



This hypothesis is confirmed in our simulations. Fig- 
ure ^ shows the distribution of time between avalanches 
for x = 10 and c = 1. It is approximately constant for 
small values and then it decays following a power law 
in more than two decades. The plateau at small values 
of At is associated to the rapid decay observed at short 
times while the tail for large At should be related to the 
long time tail. Thus, there should be some connection 
between the distribution of time between avalanches and 
the long time relaxation. 

In fact, if we assume that the relaxation is given by the 
superposition of exponential relaxations, the relaxation 
time being the time between avalanches, the temporal 
decay of N will be given by 

/>oo 

N(t) = N / dAt ra(At) exp(-t/At) . (7) 
Jo 

Then, for long times equation [7] gives the asymptotic be- 
havior 

N(t) ~ . (8) 



n(At) ~ At 



-1-/3 



(6) 



with < j3 < 1. 



Thus a = j3, the exponent of the long time tail and the 
characteristic exponent of the distribution of time be- 
tween collisions are mutually determined. This result is 
corroborated in figure ^ where it can be see that a power 
decay ~ At~ 1_Q fits quite well the distribution of time 
between avalanches for large At. Long time tails relax- 
ation and fractal distributions of time constants are the 
same phenomena manifested in different windows [|12[. 
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FIG. 5. Normalized relaxation <p(t) = N(t)/N(t = 0) of the 
number of grains in the pile for two different functional de- 
pendences p = p(N). The initial condition is N(t = 0) = N c 
in both cases. Insets show the values of a obtained from the 
fit to the power tail t~ a . 

Finally we want to mention that these simulations 
where also carried out assuming other functional depen- 
dences between the dragging probability and the num- 
ber of grains in the pile p(N). In all cases the results 
where qualitatively similar to those presented here using 
equation [| reflecting that the precise dependence is not 
determinant. For instance, in figure |^ the normalized 
relaxation for two different functional forms of p(N) is 
plotted. In both cases the long time relaxation follows 
a long time tail but with different exponents. The qual- 
itative form of the relaxation does not depends on the 
detailed form of p(N), i.e. on the detailed nature of the 
system. Nevertheless, the form of p(N) does determine 
the value of the long time tail exponent a suggesting that 
different universality classes are possible. 

We conclude that the introduction of a Bethe lattice 
representation for the avalanches and a feedback mecha- 
nism describes quite well the principal features of the re- 
laxation in sandpiles under low intensity vibrations. The 



proposed representation leads to long time tails relax- 
ation, aging, and fractal distributions of time constants, 
which are characteristic properties of the dynamics of 
complex systems out from the equilibrium. The slow 
relaxation dynamics and related properties are a conse- 
quence of the feedback mechanism, but the detailed na- 
ture of this feedback is not relevant for the qualitative 
behavior. 
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